ITERATIVE PROCESSES RELATED TO RIORDAN ARRAYS: THE 
RECIPROCATION AND THE INVERSION OF POWER SERIES 



ANA LUZON 

Abstract. We point out how Banach Fixed Point Theorem, and the Picard successive 
approximation methods induced by it, allows us to treat some mathematical methods in 
Combinatorics. In particular we get, by this way, a proof and an iterative algorithm for 
the Lagrange Inversion Formula. 



1. Introduction: beginning with a simple question. 

The results in this paper are consequences of special interpretations, as fixed point prob- 
lems, of the two classical reversion processes in the realm of formal power series: the 
reciprocation, i. e. the reversion for the Cauchy product, and the inversion, i. e. the 
reversion for the composition of series. 

The case of the reciprocation was studied in [I] and [5J. To unify our approach we also 
survey herein some of the previous results on this topic. 

The aim of this paper is to show how the (Picard) successive approximation method in- 
duced by Banach's Fixed Point Theorem (and some mild generalization) allows us to treat 
some mathematical methods in combinatorics getting so some associated algorithms. In 
particular: 

1) To construct all elements in the Riordan group as a consequence of the iterative 
process obtained to calculate the reciprocal of any power series admitting it. Presenting 
also an algorithm and one pseudo-code description for it. 

2) To construct, approximatively, the inverse of any power series (admitting it) in such 
a way that the Lagrange Inversion Formula can be first predicted and finally proved. We 
also describe the corresponding algoritm. 

For completeness we are going to recall the metric fixed point theorems we will use, see, 
for example, [1] for the first one and [10] in page 212 for the second one. 
Banach Fixed Point Theorem (BFPT). Let (X,d) be a complete metric space and 
/ : X — > X contractive. Then / has a unique fixed point xo and f n (x) — > xo for every 
x G X. 



Key words and phrases. Banach Fixed Point Theorem; ultrametric; Riordan group; Lagrange Inversion 
Formula. 
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Generalized Banach Fixed Point Theorem (GBFPT). Let (X,d) be a complete 
metric space. Suppose {/ n }neN : X — > X is a sequence of contractive maps with the 
same contraction constant a and suppose that {/„} — > f (point to point). Then / is 
a-contractive and for any point z G X the sequence {f n o ■ ■ ■ o fi(z)} > Xq, where Xq is 

n— >oo 

the unique fixed point of /. 

Our framework is the following: we consider K a field of characteristic zero and the ring 
of power series K[[x]] with coefficients in K. If g is any series given by g = J2 n>0 gn,x n , 
we recall that the order of g, u)(g), is the smallest nonnegative integer number n such that 
g n 7^ if any exist. Otherwise, that is if g = 0, we say that its order is oo. It is well-known 
that the space (K[[x]], d) is a complete ultrametric space where the distance between / and 
g is given by d(f,g) = 2 ^J_ g) , f,g 6 K[[x]]. Here we understand that ^ = 0. Moreover 
the distance between / and g is less than or equal to ^tt ; i- e - d(f, g) < i^+r, if and only if 
their n-degree Taylor polynomials are equals, T n (f) = T n (g). Finally the sum and product 
of series are continuous if we consider the corresponding product topology in K[[x]] x K[[x]]. 
See for example [7], [1] and [5] for these topics. 
In this paper N represents the set of natural numbers including 0. 
This work is motivated by the following question: 

Question 1. Can we sum the arithmetic-geometric series YlT=i kx k ~ l using the Banach 
Fixed Point Theorem?. 

We can sum the geometric series using BFPT. A visual proof of this fact can be found in 
|13J . Herein we recall an analytic proof: (the peculiar name of the following function will 
be justified later on). We consider 

h mA :(K[[x]],d) (K[[x]],d) 
t \-> xt + 1 

Since h m x is contractive, in fact d(h m ,i(ti), /i m ,i(*2)) < — cf(ti, t 2 ), we iterate at t — and 
we obtain: 
h mA (0) = 1 
h^ tl (0) = x + l 

^m,l(°) = X 2 + X + 1 

= x 3 + x 2 + x + l 

that is, 

n 

2 



As the fixed point of h mt i is the solution of xt + 1 = t, then t = and from BFPT we 
induce 



which is the unique fixed point of the function, in this case, h m> i(t) = hi(t) = xt + 1. 
Now it is natural to wonder Question [TJ 
We organize the paper in the following way: 

In Section |2] we use BFPT to a suitable function related to Question [TJ We do not answer 
the question by this way but we find an interesting arithmetical triangle. Later and using 
GBFPT we answer the question. Actually we construct the whole Pascal triangle by this 
method, see [I] and [5]. 

In Section |3] we generalize the method above, to construct the Pascal triangle, finding so 
a way to construct arithmetical triangles T(f \ g) for any pair of series / and g with non 
null independent terms. Using the usual product of matrices, we identify the well-known 
Riordan group, see [I]. 

In the procedure described above we obtain a new parametrization of the elements in 
the Riordan group and so a new notation different from the usual ones. In Section 4 we 
try to justify the use of our notation alternatively to the usual notation. Our method of 
construction and our notation allow us to explain easily a way to add and delete columns 
suitably in a Riordan array to get another one. For a concrete kind of triangles, those 
denoted by T(l | a + bx), we can calculate the inverse only adding adequately new columns 
to those triangles. In fact it is an elementary operations method. We end this section giving 
expressions for the so called A and Z sequences of a Riordan array. These expressions are 
given in terms of our notation and related to the inversion in the Riordan group. 
In Section we give the main new results of this paper. We display an algorithm to 
construct the inverse of a series and we show the relation with the Lagrange Inversion 
Formula. In particular we prove that Banach Fixed Point theorem gives rise to the Lagrange 
Inversion Formula. 

2. TWO ANSWERS: A CURIOUS TRIANGLE AND AN ITERATIVE METHOD 

To answer the previous question, in this section we are going to recall lightly some examples 
and tools widely studied in [I] and [5]. 

It can be easily shown, see page 2257 in [5], that there are not a first-degree polyno- 
mial with coefficients in K[[se]], f(t) = g(x)t + h(x), and any point, Xq, such that its 
n + 1-iteration coincides with the partial sum of the arithmetic-geometric series, that is: 




k=0 



$f,x / Vn, r+ Vo) = £Lo(* + 1 > k - 



3 



In view of this, we are going to iterate a polynomial whose fixed point is the sum of the 



(l-x)' 2 



Since the equality t 



arithmetic-geometric series, that is Y2t=i kx k ~ 
be converted to t — 1 + (2x — x 2 )t, we consider the polynomial f(t) 
coefficients in and we initiate the iteration process at t — 0: 

/«') = 1 

2x-x 2 

2x + 3x 2 -4x 3 + x 4 
2x + 3x 2 + 4x 3 -llx 4 + 6x 5 - x 6 
2x + 3x 2 + Ax 3 + 5x 4 -26x 5 + 23x 6 



l + (2x- x 2 )t, with 



/ 2 (0) 
/ 3 (0) 
/ 4 (0) 
/ 5 (0) 
/ 6 (0) 



8x 7 



,8 



1 + 2x + 3x 2 + Ax 3 + 5x 4 + 6x 5 -57x 6 + 72x 7 



x 
39x 8 



10x 9 



,10 



We know that the sequence of iterations converges to the sum of the arithmetic-geometric 
series, because / is contractive in (K[[sc]], d), but we can observe that in each iteration the 
partial sum of such series appears plus a remainder. We want to control the difference with 
the partial sum. To do this, we display the coefficients of the remainder as a matrix, that 
is: 

\ 



( -1 






-4 


1 




-11 


6 


-1 


-26 


23 


-8 


-57 


72 


-39 


-120 


201 


-150 



1 

10 -1 
59 -12 1 

V : : : 
Observing this matrix we recognize: 

(1) The rule of construction is similar to that of Pascal triangle: each element is twice the 
above element minus the element above to the left side, that is: a n j. = 2a n _i £ — a„,_i 

(2) The elements in the first column are Eulerian numbers except for the sign. 

(3) The sum of the elements in any row are triangular numbers with negative sign. 

(4) For every element, the sum of all elements in its row to the right and all elements above 
in its column is zero. That is, Y^k=i a kj + Ylk=j+i a ^ = 0- 



(5) The general term is a n j 



n 



+j-i+Eti(-i) fc (tS^)2" + - 



-2k 



Etc. 



For an exhaustive development of this triangle see [5]. 

The above approach, using BFPT, does not give us an exact answer to our question. To 
find an adequate answer we consider the GBFPT. This is our way to do this: 
For computability facts we consider the sequence of functions, with polynomial coefficients, 
given by ho^t) = xt, h 12 (t) = xt + x, h 2j 2(t) = xt + x + x 2 , h 3t2 (t) = xt + x + x 2 + x 3 , 



m—l 



h m ,2{t) = xt + x^2 



X 



k=0 



Each function h m , 2 is —contractive, so {h mi2 } is an equi-contractive sequence of one-degree 
polynomials that converges to h 2 (t) — xt + y^, i. e. : 

x 

{h m ,2} — > h 2 (t) =xt + —— 

-L Ob 



It is easy to see that the crossed iterations induced by GBFPT are just the corresponding 
partial sums of the arithmetic-geometric series: 

M°) = 

^1,2(^0,2(0)) = x 

^3,2(/i2, 2 (/il,2(/iO,2(0)))) =X + 2X 2 + 3X 3 

Now using again GBFPT we obtain that, since xt + = t =>- 



(l-x): 



-, then 



(h 



m,2 



V)(0) 



.7; 



1 - x 



these crossed iterations at zero converge to the unique fixed point of h 2 , that is, the sum 
of the arithmetic-geometric series. So the answer to our question is yes if we are allowed 
to use the generalized version of the BFPT. 
Recall that the Pascal triangle is given by: 



1 

3 1 

6 4 

10 10 

15 20 



1 

5 

15 



(0) (?) (3) (3) (2) (5) (e) 



CO 



(1-x) 2 (1-x) 3 (1-z) 4 (l-a:) 5 (l-x)b (l-x)7 



(l-x) n 



We have just constructed the first two columns of Pascal tringle using BFPT. In fact we 
needed only BFPT to construct the first one and GBFPT to get the second one. The 
main observation is that we can follow this iterative procedure to construct all columns. 
For example we can repeat the process to construct the third column. To get this goal, we 
interpret the above equicontractive sequence h m ^ in the following way 



m— 1 



h m , 2 {t) = xt + x^2x k = xt + xT m _ x ,\ 



fc=0 



where T m _\^ is the m — 1 degree Taylor polynomial of the first column (which is the 
geometric series). So in a similar way we consider the following equicontractive sequence: 

m— 1 

h m ,3(t) = xt + x^2kx k = xt + xT m _ lj2 



fc=0 



where T m _i ; 2 is the Taylor polynomial of the second column (which is the arithmetic- 
geometric series). So, as one can easily prove, the crossed iterations for this sequence 
coincide with the partial sums of the third column: 

Mo) = o 
^1,3(^0,3(0)) = 0, 

MMM°))) =x 2 , 
MMMM0))))=z 2 + 3rr 3 , 

MMMMMo))))) = x 2 + 3x 3 + 6x\ 

Using once more GBFPT, we obtain that these crossed iterations converge to the unique 
fixed point of the limit function h 3 (t) = xt + Xjj3^2- Since 



2 

x x 



h 3 (t)=xt + x— — =>- t 



[l-xf iX-xf 
then 



(V3---M(o) = £ ( W 

1 — n V / 



x 2 



X) 3 



k=0 

Actually, as we said before, we can construct every column of Pascal triangle using this 
process: 

Proposition 2. For n > 2, the n-column in Pascal's triangle is obtained from the (n- 
l)-column applying the crossed iterations in GBFPT to the sequence {h rntn } m ^m where 

hm,nif) Xt + xT m —\ n —i 

being T m „ 1 „_ 1 the (m-1) -Taylor polynomial of the (n-l)-column. 



3. THE GROUP OF ALL ARITHMETICAL TRIANGLES T(f | g). 

Now we generalize the previous iterative method for any pair of series / = Xln>o Z" 3 -™ an d 
g = J2n>o9 nXn sucn that /o ^ and g 7^ 0. We construct the following arithmetical 
triangle T(f \ g), in this notation the Pascal triangle is T(l | 1 — x), where the role of series 
1 is played by the series / and the role of 1 — x by g. 



fo 








fl 


do,o 






/ 2 


di,o 


do,i 




fs 


^2,0 


d 2 ,i 


^2,2 


f 


L 

9 


xf 
9 2 

6 


x 2 f 
9 3 



In [1] we interpreted the calculation of — as a fixed point problem. Consider the sequence 

h m ,i{t) = T m (— — - J t + T m (— 

where T m (f) is the m degree Taylor polynomial of /. Observe that the sequence of crossed 
iterations has as its limit the unique fixed point of hAt) = — — —t + — , that is: — . It is 

go go g 

the first column of T(f \ g). 

To construct the second column and the next ones we consider the equicontractive sequences 

h (+\-t( 9o~9 \ , T ( xU ~ 2 f 

V go J Kgog 11 ' 1 

their corresponding limits are h n (t) = [j^-J t + x ^ g ogn -i J whose corresponding unique 
fixed points are 

t = L 

g n 

The series t n is just the n-column of our T(f \ g). 

Theorem 3. Let f = f n x n and g = g n x n with g ^ 0, then the Riordan matrix 

n>0 n>0 

T(f | g) = (d ntk ) is given by 
if k = 0, then 

, _ fo , 9l, <?2 , 9n , . fn 

"0,0 — , "n,0 — "n-1,0 "n-2,0 ' ' ' "0,0 H 

go go go go go 



if k > 0, then 



j #1 , <?2 , g n -k , G?n-l,fc-l 
a n,k — "n-l,fe "n-2,fc ' ' ' "fc.fc + 



go go go go 

This theorem gives us the following algorithm to construct for columns any Riordan matrix: 
Algorithm 4. Given / = f n x n and g = '^^g n % n with g ^ 0: 

n>0 n>0 

Step 1: Calculate the first column d n . 

i _ fo , gi, gi , 9nj.fn 
"0,0 — — > «n,0 — "n-1,0 "n-2,0 ' ' ' "0.0 i 

go go go go ' go 

Step k: Calculate the fc-column, d n> k, using k — 1-column. 

, 9l , g2 , 9n-k j dn-l,k-l 
«n,fc — "n-l.fc "n-2,fc " " " "fc.fc H 

go go go go 



We can write this algorithm in an informal pseudo-code: 
READ (f,g,n) 
SET (d,aux) 

CAL CULA TE d[0, 0j=f[0j/g[0j 
% We calculate the first column 
FOR i=l to n 
FOR k=l to n 

CALCULATE aux(k,i)=g[i-k]*d[k,0] 
END 

CALCULATE d(i,0)=l/g[0]*(f[i]-SUM(aux(:,i)) ) 
END 

% We calculate the remaining columns 
FOR j=l to n 
FOR i=l to n 
FOR k=l to i 

CALCULATE aux(k,i)=g[i-k]*d[k,j] 
END 

CALCULATE d(i,j)=l/g[Q]*(d(i-l,j-l )-SUM(aux(:,i)) 

END 

END 

PRINT(f,g,d) 

So to construct the arithmetical triangle T(f \ g) it is enough to know the ordered pair of 

series / and g, i. e. the data, and the algorithm of dividing two series. Every column is 

/ / 
constructed by the same rule as that in — but the coefficients of — are replaced with the 

9 9 

coefficients of the previous column. Except for the first column, here we need and auxiliary 

column, the coefficients of /. 

We can consider the matrix T(f \ g), like in Linear Algebra, as the associated matrix to a 
K-linear continuous function, see [3]: 

T(f\g): (K[[x]],d) (K[[x]],d) 

h h. T(f\g){h) = £h(z) 

Using the classical definition of composition of maps and the behavior of the associated 
matrix, we can easily find the formulas for the product and the inverse for these triangles. 



T{h I 9i)T(f 2 | 92) = T M 



9\9i 



(T(f I = T 



/(w-i) g(co-i)J> 



uo o Id 



0J O Id 



X 



So if we consider the set of the all arithmetical triangles with /o ^ and go 7^ and 
the usual product of matrices we obtain a group. Actually this group is the well-known 
Riordan group. 



4. On the T(f \ g) notation. 

We have received some critics about our notation. Someone could think that our notation 
is, in some sense, cumbersome. Of course it depends strongly on the way you approach or 
you run into this group. In this section we are going to give some reasons why our notation 
could be very adequate. The basic formula relating our to the classical notation is 

70) % 

' 1 

9, 



h 



g{x) 



(d(x), h(x)) 

The fundamental equality with our notation is: 

(1) T(f\g) = T(f\l)T(l\g) 

This equality in the other notation is: 

77 



T(f I 9) = (d 



i,j)i,j>0 



X 



X 



= (f(x),x) 
g{x) g{x)j \g{x) g{x) 

In ([1]) we can see that every element of the Riordan group [9] can be expressed by means 
of the product of a lower triangular Toepliz matrix whose columns are the coefficients of 
series /, shifted conveniently, the matrix T(f | 1), and a renewal array, the matrix T(l | g) 
described by Rogers in [5] . These last kind of matrices are really similar to the Jabotinsky 
matrices, see [3] . We want to point out that the structure of every element of the Riordan 
group is essentially in the structure of the matrix T(l | g). For example, to know a closed 
formula for the general term of T(l | g) gives us at once a closed formula for the general 
term in T(f \ g). This matrix T(l | g), for us, is intrinsically related to the calculation of 

— , which is its first column. 

g 

A comparative table of both notations is given: 



Name 


(d(t),th(t)) 


T(f|g) 


Identity 


(1,*) 


T(l | 1) 


Pascal 


( 1 , M 
\i-t i-t j 


T(l | 


Appel subgroup element 


(d(t), t ) 


T(d | 1) 


Associated subgroup element 


(i,th( t )) 


T {h 7) 


Bell subgroup element 


(d(t),td( t )) 





A curious and symbolically important, for us, property of our notation is the way to give, 
by means of the parameters, the natural powers of the Pascal triangle: 



Proposition 5. For every n GM, we have T n (l | 1 — x) — T(l | 1 — nx) 



Proof. Let us proceed by induction. 

For n = 2: T 2 (l | 1 - x) = T(l | 1 - x)T(l | 1 - x), u 



So 

1— x 



T 2 fl I 1 



x 



T(l|(l-x)(l 



.r 



T(l | 1 - 2x) 



I — x' 

Suppose that T n ~ l {\ | 1 — x) — T(l | 1 — (n — then 

T n (l | 1 - x) = T(l | 1 - x)T n ~ l (l | 1 - x) = T(l | 1 - z)T(l | 1 - (n - l)z) 



T(l | (1 - x)(l - (n - 1) )) = T(l | 1 - na;) 



x 



□ 



Another reason why for us our notation is natural, is related to the way we begun to study 
these topics. One of the first things we did was to find our curious triangle described in 
Section [2j From our notation, the description as a Riordan array is: 

\ 



/ -1 










-4 


1 








-11 


6 


-1 






-26 


23 


-8 


1 




-57 


72 


-39 


10 


-1 


-120 


201 


-150 


59 


-12 



T 



(1 



2x - 1 



x) 



This notation resembles both the problem we were treating and the algorithm of construc- 
tion. 

Another thing we can describe easily with our notation is the fact that, with our construc- 
tion method by columns, we can add new columns to the left for every element of the 
Riordan group to obtain again a Riordan array intrinsically related to the initial one, for 
example: 



(\ 










2 


-1 








3 


-4 


1 






4 


-11 


6 


-1 




5 


-26 


23 


-8 


1 


6 


-57 


72 


-39 


10 -1 


7 


-120 


201 


-150 


59 -12 1 


V 








■ '■J 



2x 



1-x) 



2x 



Note that we added a new column to the left and look at the way the parameters changed 
in our notation. 
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In general we can construct a family of new Riordan matrices closely related to it. For 
example by definition of Riordan array we get 



T(fg | g) 



( fo 

fi d fi 

f 3 rf 2 ,0 ^2,1 ^2,2 

/ 4 rf 3i0 d 3j i rf 3i2 ^3,3 

V ; ; ; : ; 



/ (-15 

.9 



( di,i 

dl,\ C?2,2 

^3,1 ^3,2 ^3,3 

d A i d^2 ^4,3 ^4,4 

^5,1 ^5,2 ^5,3 ^5,4 G?5,5 

V ; ; ; ; ; 



\ 



where / = ^ / n x n and T(f \ g) = (d n ,fc)n,fceN- 

n>0 

By the same way we observe that T | g j for m G N is the matrix obtained from T(f \ g) 
by deleting the first m-rows and m-columns. Moreover T(fg m \ g) is the unique Riordan 
matrix with the property that by deleting the first m-rows and m-columns from it, we 
obtain T(f \ g). In fact it can be easily proved the following: 



Proposition 6. Let T(f \ g) = (d n ^) n ,kem be a Riordan matrix, and m e Z then 
T(fg m | g) = (d n ,k)n,keN, with d n>k = [x n - k ]fg m - k - 1 . Where [x j ]S stands for the j- 
coefficient of the formal power series S . 



We can always embed any T(f\g) in a bi-infinite lower triangular matrix. 
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For our example we have 



/ 
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-23 


-6 


-1 























26 


11 
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1 




















-5 


-4 


-3 


-2 


-1 

















-4 


-3 


-2 


-1 





1 














-3 


-2 


-1 





1 


2 


-1 











-2 


-1 





1 


2 


3 


-4 


1 








-1 





1 


2 


3 


4 


-11 


6 


-1 








1 


2 


3 


4 


5 


-26 


23 


-8 


1 



This construction is very nice in the case of Riordan matrices of the kind T(l | a + bx), 
(treated as change of variables in [5]). For these matrices we have 



/ 



\ 



i (_s_\ 4 

a \l — bxJ 

a VI— bx J 

- ( — ? 

a \l—bx/ 



a 



1 

a l—bx 
1 
a 

1 n-bx \ 
a \ a J 

l n-bx \ 2 

a \ a J 





a 3 




















<- 


3a 2 b 


a 2 



















3ab 2 


2ab 


a 














<- 


b 3 


b 2 


b 


1 











4- 














1 

a 








4- 














b 

a 2 


1 

a 2 





4- 














b 2 
a 3 


2b 

a3 


1 

a 3 




1 


1 


1 


I 




1 






(a + bx) 3 


(a + bx) 2 


(a + for) 


1 


1 


1 


1 




a+bx 


(a+fcr) 2 





V 



Note that, as in the Pascal triangle, we can see the above matrix in the following way 

\T- 1 (l\a + bx)+ 

T(l\a + bx) 

Where ^T _1 (l | a + bx)t is ^T _1 (l | a + bx) placed in the same way as the inverse 
of the Pascal triangle is placed in the Hexagon of Pascal in page 194 in [2]. Note that 
T _1 (l | a + bx) = T(l | ^-—^)- This gives us a method to calculate T _1 (l | a + bx) by means 
of elementary operations. 
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Our algorithm of construction do not need the so called A and Z sequences of a Riordan 
array, see [B] and [UJ. But, in our notation, they appear in the expression of the inverse of 
the Riordan array T(f \ g) giving us specially aesthetic formula: 

Proposition 7. Let f = f n x n and g = g n x n be two formal power series with 

n>0 n>0 

/o ^ and go ^ 0. Suppose that A and Z represent the A-sequence and the Z -sequence, 
respectively, ofT(f \ g). Then 



(i) T-\l\g)=T{l\A) 
(") T-\f\g)=T[f Q {A-xZ)\A 



Proof, (i) From Theorem 1.3 of [TT], the A-sequence is the unique series, with A(0) ^ 0, 
such that — = A ( — Y So — = xA ( — ) . If co = — then co = xA (co). On the other hand 

g \gj 9 \gj 9 

as co' 1 o to = x and to = — then x = cog, composing with co~ l we get co' 1 = xg(co~ l ) and 

9 

g^uo' 1 ). So composing with co~ l but now in co = xA (co) we get x = co A(x) then 



CO' 1 



X 

,-1 



1 = —A(x) so 1 = g(to- 1 )A(x) then A(x) = — Since T'^l | g) = T(l | 

x g[co i ) 9( ' 

T(l | A). 

(ii) From Theorem 2.3 in [6] we obtain that the Z is determined by the equality 
io~ x Z = 1 — TT - . From here we get —. — -- = ^j-(A — xZ). So 

□ 

Corollary 8. 



T- 1 (/|( 7 ) = T(1|A)T 



5. Lagrange inversion formula via Banach fixed point Theorem 

In the previous section we showed that to calculate the inverse of T(l | g) we need, in 
particular, to calculate co' 1 , where co = ~ and then co' 1 = xgiuo' 1 ). So we consider the 
function F : iK[[x]] — )■ xK[[i]] defined by F(y) = xg(y). Here xK[[x]] represents the series 
with null independent term. This function is —contractive since 

d(F{ yi ),F(y 2 )) = 2w[xg{y l_ xg{y2) ) = 2 -( 9 (,i)- 9 fe))+i ^ 2 d ( yi,V2) 

13 



The domain, of F is the closed ball, in (K[[x]], d), whose center is the series 5 = 

and the ratio is |. Consequently our domain is also complete with the relative metric. So 

the unique fixed point of F is u^ 1 = (^j and BFPT can be applied. 
The BFPT gives us a theoretical iterative process to calculate u^ 1 . To convert this method 
into an effective approximation process we first note that the relation d(Si,S 2 ) < ^ts+t 
means that the m degree Taylor polynomials of both series are equals, that is T m (Si) = 
T m {S 2 ). Then we obtain the following algorithm: 

Suppose g = ^2 n>0 g n x n - We begin to iterate at S — 0. F(0) = g x. This means that 
T^u; -1 ) = g x. Using the |-contractivity of F we get T 2 (F(g x)) = T 2 (uj~ 1 ). Since 
F(g x) = g x + gogix 2 H — • we obtain T 2 (u;~ 1 ) = g$x + g^gix 2 . Similar arguments allow us 
to prove that T 3 (F(g x + gogix 2 )) = T 3 (o;~ 1 ). The above construction can be summarized 
in the following (the notation is as above): 

Proposition 9. 

T m {uj- 1 ) = T m (F(T m _ 1 (F(. • • (FCT^O)))) • • • )))) 

Following this process we get 

Ti(u- 1 ) = g x 

T 2 {uj- 1 ) = g x + g gix 2 

T 3 (u~ r ) = g x + g gix 2 + (g g 2 + g 2 g 2 )x 3 

T 4 (w~ 1 ) = g x + g g ± x 2 + (g gj + g$g 2 )x 3 + (g gf + 3g 2 g x g 2 + glg 3 )x A 

T 5 (w~ 1 ) = g x + g gix 2 + (g g 2 + glg 2 )x 3 + (g gf + 3g 2 g t g 2 + glg 3 )x A + (g gf + 6g 2 g 2 g 2 + 

Zgfal + 4g$gig 3 + g^)x 5 

If we recall the Cauchy powers of the series g: 

9(x) = go + 9ix + g 2 x 2 + g 3 x 3 + g^x 4 H 

g 2 (x) = g 2 Q + 2g gix + (2g Q g 2 + g 2 )x 2 + (2g g 3 + 2g x g 2 )x 3 ■ ■ ■ 

g *(x) = gl + Sgfaix + 3(g g? + gog 2 )^ 2 + (6# #i#2 + 3g$g 3 + g\)x 3 + ■■■ 

g\%) = go + ^gfoix + (4^ 2 + %o#i> 2 + 4 (gog? + 3gggig 2 + gogs)^ 3 + • • • 

g\x) = gl+hg 4 g 1 x+(hg 4 g 2 +l^g 3 g 2 )x 2 +{2Qg 3 g l g 2 +\{}g 2 ^^ + 6g 2 g?g 2 + 2g ! gl 

+4gogigs + gog4)a: 4 H 

comparing adequately the coefficients of u^ 1 and the powers of g we obtain the next 



relationships: 


[x]u- 1 -- 


= [x°]g 


[x 2 ]u- 1 


= \W 


[x 3 ]uj- 1 


= \w 


[x^co- 1 


= \w 




= \W 
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These equalities allow us to predict and motivate the classical Lagrange Inversion Formula, 
see [T2] page 36: 

[x n+l ]uj- 1 = — '— [x n ]g n+l , with ui = - 
n + 1 g 

From now on we denote Tj = Tj(u~ l ). To show how this process works note that 
F{T n ) =x{g + 9l T n + g 2 Tl + • ■ ■ + <? n T* + ■ ■ ■ ) = 
= T n + ( 9l [x n }T n + g 2 [x n ]T 2 n + ■■■+ g n [x n ]T^)x n+1 + S n+2 with S n+2 G x n+2 K[[x]] 

So 

n 
k=l 

Because [x n ](c<;~ 1 ) A: = [x n ](T n ) k for any k < n. Suppose now that we know 

n[x n ]{uj- l ) k = k[x n - k ]g n for k < n 

then 

1 n 1 1 

[x n+l }F(T n ) = [x n+1 ]uj- 1 = -J2 k[x k ]g[x n - k }g n = -[x n - l ]g'g n = — [x n ]g n+l 

n k=1 n n + 

Note that in the above development we need to know n[x n ](u~ 1 ) k = k[x n ~ k ]g n for k < n. 
In fact we can give a proof of all above using essentially the fact that u" 1 is a fixed point 
of certain contractive function. 

Theorem 10. (Lagrange inversion via Banach Fixed Point Theorem) Let IK be a field 
of characteristic zero. Suppose that u is a formal power series in K[[x]] with w(0) = and 
u/(0) ^ 0. Then 

n[x n ](u^) k = k[x n - k ] [-) for n,keN 

x 

Proof. Let g = — . So [x°]g ^ 0. As proved before a; -1 is the unique fixed point of the 

CO 

^-contractive function F : xK[[x]] — > xK[[x}} defined by F(y) = xg(y). Iterating at y = 
we get 

[x 1 ]^- 1 = [x l ]F(0) = [x°]g 
If k > 1, note that [x 1 ](o;~ 1 ) fc = and [x^^g = and then 

[x 1 ]^- 1 = k[x l - k ]g 

Let proceed by induction on n. Suppose that 

j[x j ](u' 1 ) k = k[x j ~ k ]g j for j < n, k > 1 

Note that there are actually only a finite number of suppositions on k. Because if j < k, 
then [x^^gi = = [^^(a;^ 1 )^. Then the equality holds trivially. 
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Since u 1 = xg(uu x ), then for any k (uu v ) k = x k g k (uj r ). Consequently 

n+l-fe 

[x n+1 }{uu- 1 ) k = [x n+1 ]x k g k {u>- 1 ) = [x n+1 -V(^ _1 ) = Yl \Ag k [x n+l - k ]{uj- l y 

3=0 

by the induction hypothesis 

n+l—k 

[x n + i ]{uj -i )k = J- 3lx>]g k [x n+1 - k -i]g n+1 - k 

Let us call h = g k 

n+l—k n+l—k 

k" +1 ](^ 1 ) fc = n + 1 _ fc J[x j nx n+1 - k ~ ] ]h n± ^ = n + 1 _ k £ [xi- l \h'[x n ^- k -i\h 

j=0 j=l 
n— k 

= — - — -VHfe'b" 4 -^ 5 ^ = — - — -ix^m'h^) = 

n+l-k ^ l 11 1 n+ 1 - k l JV 7 

j"=0 



= l -^[x n - k ] ( — h^) = J^[ X ^- k ] h ^ = J^[ X ^- k ]g^ 

n+l-k l J V^ + 1 / n + V J n + l L J 

□ 

The development above gives us the following algorithm to calculate the coefficients of the 

x 

compositional inverse of u = —: 

9 

Algorithm 11. Given g = ^^gjx^, with g ^ 0. Given F(y) = xg(y), y G a;K[[a;]]. 

j>o 

step 1: (Initial.) T x = g x. 

step k (2 to n): Calculate the Taylor polynomial of order k of F(T fe _i) 

We can write this algorithm in an informal pseudo-code: 
READ (g,F,n) 
SET T 

CALCULATE for i from 2 to n do 

T[i]:= convert ( series (F( T[i-1], x=, i ) ),polynom ); 

end 

PRINT T 
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